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Preface 

This  is  the  first  in  a  series  of  notes  on  various 
topics  in  magneto-hydrodynamics,   A  large  part  of  this 
material  has  been  presented  in  informal  lectures  at 
irregular  intervals  since  1951;.   A  set  of  four  mimeo- 
graphed notes,  representing  some  of  the  early  develop- 
ments, was  issued  in  July  195^  with  the  purpose  of 
acquainting  a  small  group  at  NYTJ  with  certain  aspects 
of  the  subject.   Meanwhile,  a  certain  amount  of  general 
ixaterest  has  been  exhibited  in  these  notes,  and  it  was 
felt  to  be  worthwhile  to  rewrite  the  original  notes  in 
a  somewhat  more  readable  form,  at  the  same  time  extend- 
ing and  elaborating  them  greatly^   The  original  four 
notes  are  contained  in  numbers  1,  3,  ^4-,  7  of  this  series. 
Although  many  of  these  reports  are  unconnected,  an  attempt 
will  be  made  to  supply  enough  background  material  to  make 
the  series  to  some  degree  self-contained. 

The  status  of  the  subject  of  magneto-hydrodynamics 
is  somewhat  paradoxical.   It  involves  in  combination, 
various  aspects  of  fluid  dynamics,  electromagnetic  theory, 
kinetic  theory,  and  occasionally  relativity.   Most  of 
this  was  available  to  Maxwell  and  he  had  the  material  to 
develop  the  subject  to  a  considerable  extent,  but  only 
recently  has  it  been  attacked  significantly.   Another, 
and  probably  related,  consideration  is  that  the  companion 
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subjects  of  fluid  dynamics  and  electromagnet! sm  have  a 
well  developed  experimental  backgroiond,  while  relatively 
little  is  known  experimentally,  in  magneto-hydrod3majnics. 
Consequently,  the  development  of  magneto -hydro dynamic 
theory  is  rather  abstract,  principally  devoted  to  the 
solution  of  rather  idealized  special  problems,  and  the 
directions  of  development  are  guided  to  a  large  extent 
by  what  has  been  done  in  fluid  dynamics  and  electromag- 
netics.  There  do  exist  important  problems  in  astrophysics, 
in  properties  of  the  earth's  atmosphere,  and  in  gaseous 
discharges,  but  these  are  usually  so  difficult  as  to  be 
unapproachable  in  the  framework  of  a  mathematical  theory. 
Controlled  experiments,  with  only  a  few  effects  present 
at  a  time,  are  very  scarce. 

This  first  note  is  largely  a  matter  of  bookkeeping 
and  can  be  motivated  only  by  pointing  to  the  various 
applications  made  later  on  or  by  reference  to  kinetic 
theory.   Various  forms  of  fluid  conservation  equations 
are  listed  for  future  reference.   They  are  presented  in 
forms  too  general  to  be  of  specific  mathematical  use,  but 
they  do  serve  the  useful  purpose  of  correlating  and  categor- 
izing the  many  specific  fluid  descriptions  which  will  be 
used. 

This  note  was  prepared  with  the  assistance  of 
Ignace  Kolodner  and  Albert  Blank. 
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General  Fluid  Equations 

Any  macroscopic  picture  of  a  fluid,  electrically 
conducting  or  not,  is  based  on  conservation  laws,  usually 
of  mass,  raoinentum,  and  energy,  and  occasionally  of  others. 
It  is  Important  to  realize  several  things.   First,  these 
are  "lack  of  conservation"  more  than  "conservation"  equa- 
tions.  For  example,  the  energy  of  an  element  of  fluid 
(or  of  the  whole  fluid)  is  usually  not  a  constant,  but  is 
assumed  to  vary  in  certain  definite  ways.   Second,  the 
form  taken  by  these  equations  is,  macroscoplcally  speak- 
ing, pure  postulation.  One  is  governed  by  certain  general 
considerations  of  what  has  been  found  to  be  useful,  but 
no  question  of  derivation  can  arise  iinless  referred  to  a 
more  primitive  basis  such  as  motion  of  particles.  This 
remark  applied  also  to  the  t he rmo -dynamic  formalism  which 
is  microscopically  derivable  but  macroscoplcally  postu- 
lated.  Finally,  in  a  strictly  formal  sense,  many  of  the 
terms  included  below  are  redundant.   For  example,  in  the 
momentum  equations,  we  find  terms  X  and  ^P  ^/ix^% 
Clearly,  it  would  be  completely  general  to  omit  ^P  "^/3x"' 
and  claim  that  the  force  term  X   sometimes  takes  the  form 
of  space  derivatives  of  auxiliary  quantities,  and,  con- 
versely, it  is  easy  to  see  that  a  force  term  X  can  always 
be  put  into  the  form  ^P  'V^x  by  suitable  choice  of  P  •'. 
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This  separation  into  terms  is  significant  only  because 
the  conservation  equations  are  Incomplete,  mathematically 
speaking,  and  the  supplementary  relations  which  will 
eventually  have  to  be  introduced  will  be  of  several  kinds 
for  which  the  various  forms  such  as  ^p■'■V^^  and  X"*"  are 
more  convenient.   It  may  be  remarked  that,  even  from  a 
molecular  standpoint,  the  separation  into  JP  ^/d7i,   and 
X  turns  out  to  be  somewhat  arbitrary.   Insofar  as  the 
equations  listed  below  are  sufficiently  general  to  cover 
most  cases  of  normal  occurrence,  they  cannot  be  complete 
but  must  be  supplemented  by  relations  which  specify  the 
particular  nature  of  the  fluid  under  study;  these  rela- 
tions will  range  from  functional  relations  (e.g.,  an 
equation  of  state)  to  differential  equations  (e.g.,  rela- 
ting stresses  to  velocity  gradients,  or  relating  a  force 
term  to  a  magnetic  field  vector  which  in  turn  is  coupled 
to  other  of  the  fluid  variables  through  a  set  of  auxili- 
ary electromagnetic  differential  equations). 

By  the  inclusion  of  general  terms  like  X  ,  we  obtain 
complete  and  universal  validity  for  these  equations  (e.g., 
by  defining  X   as  the  force  required  to  make  the  momentum 
equation  correctl),  but,  by  the  same  token,  they  cannot 
serve  more  than  a  bookkeeping  piJrpose.   This  use  is 
Important,  however,  because  of  the  great  variety  of  fluid 
models  which  may  be  used  to  describe  a  conducting  medium. 
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1»   Notation 

We  consider  a  fluid  consisting  of  n  homogeneous 

components,  each  with  its  ovm  Intrinsic  properties  (such 

as  molecular  mass,  charge,  etc.)  and  each  with  its  own 

thermodynamic  state  variables. 

Molecular  properties: 

m^  =  mass  ,    . 

r 

e  =»  charge  (signed)        r  =  l,  ••»,  n 

'r    r  r 
Molecular  structure  is  referred  to  only  to  establish  the 

connection  between  mass  density  and  charge  density  below. 

State  variables: 

V   =  number  density 

p   =  m  V   =  mass  density 

q   =  e  V   =  charge  density  (signed) 

q  /p   =  Y  -   constant  for  a  given  fluid 

u  =  macroscopic  flow  velocity 

p  u  =  momenttim  per  vmit  volume 
^r  r  ^ 

e   =  Internal  energy  per  unit  volume 

12 
e   =  total  energy  per  imit  volume  ~  ®v,  "*"  p  P-n^-r 

T   =  temperature 

p  =  pressure 

^i,  =  chemical  potential 

Supplementary  variables: 

P^''  =  stress  tensor,  pj''  =  Pj''  -  p^..^^*^ 
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Q^  =  heat  flow 

X  =  force  per  unit  volume 
r         ^ 

C   =  energy  source  per  unit  voliitne 

(T  =  mass  source  (chemical  creation)  per  unit  volume 
r 

Cw=   momentum  source  (chemical)  per  unit  volume 
r  r 

Internal  energy,  stresses,  etc,  are  defined  to  be  those 
observed  moving  with  the  individual  fluid,  i.e.,  at  the 
velocity  u  ♦ 
2«   Conservation  Equations 

The  follovdng  "conservation"  equations  are  postulated:-"- 

-W  "-^^    (Pr"r)  =  ^r  ^^^"^^ 

^^^  {   it    ^Pr^r^  "■  O  ^^v'V^i  "-   ^r^^  =  ^J  ""  ^r^r   ("'Omentum) 


ae 

^      ax 


""  -"  ^   (^r^r  ""  ^r^r'^'  +  Q^)  =  ^  (energy) 


Using  V  =  u  -  u,  Jl  Pt.%  =  0  we  define  the  following 
gross  (single-fluid)  variables: 
V  =  I  v^ 

p  =  Z  pj, 
q  =  I  Op 

pu  =  I  Pr^r 


-;>■  The  summation  convention  is  used  for  repeated  tensorial 
indices,  but  not  for  indices  distinguishing  the  fluids. 
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e  =  I  (e^  +  1/2  p  u^)  -  1/2  pu^  =  I  (e^  +  l/2  p  v^ ) 


e  =  51  e   =s  e  +  1/2  pu^ 


r 

X  =  I  X^ 


r 
-r 


The  V  are  the  relative  or  diffusion  velocities.   The 
r 

single  fluid  internal  energy,  stresses,  etc.,  are  defined 

with  respect  to  the  mass  flow  velocity,  u,  and  are  not 

merely  sums  of  the  corresponding  quantities  of  the 

constituent  fluids. 

Prom  the  conservation  of  mass  and  an  appropriate 

definition  of  CJ  , 

r 

I  (y^(^  =  0.  ;;:•:,  , 

Adding  the  individual  conservation  equations:     _.  ,, 

^  +-^  (pu^)  =  0  ,.  ,(  • 


ax- 


(2) 


4r  (pu^)  +  -K   (pu-^uJ  +  ?^h   =  X^ 
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It  is  sometimes  convenient  to  consider  the  following 
(Lagrangian)  form  of  the  momentiain  equations: 


(3) 


ih) 


— -  +  w^   — ? 


)t 


^x 


dp 


ij 


^:c 


J 


=  x^  +  (T   (coi  -  ui) 


and  also   the  difference  momentvun  equations 


^v 


^u 


>P. 


^  ^  ^  ^        Pr  )x 


1.1 


1  ap 


ij 


p  JJ 


X-^        ^i        (5-  .  . 

=  -^  -2L  +  -I  (u,=^  -  u^) 

Pp        P  Pp  I'  r 

Equations  (3)  and  (U)  are  obtained  by  combining  the 
momentiAm  and  mass  equations  of  (1)  and  (2)  respectively « 

The  conservation  equations  (2)  are  equivalent 
(assuming  differentiability  of  all  quantities)  to  the 
integral  relations 

1^  J  pdV  +  i  pu^dS  .  =  0 
D       S 


(6) 


1^  j  pu^dV  +  6  (pu^u^  +  P^'^)dS,  =  J  X^dV 


1^  J  edV  +  a  (eu^  +  u^P^^  +  Q^  )dS  =  /(fdV 


edV  + 
D       S 
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The  surface  integrals  are  taken  over  the  boundary, 

S,  of  the  fixed  domain,  D.  The  Lagrangian  form  of  (6) 

for  a  domain  D(t)  which  is  carried  along  by  the  fluid 
at  the  velocity  u  is  , 

4t  /  pdv  =  o 

D(t) 


it    J   p-'^^  V     ^'^^sj  =  J    ^'^ 

D(t)        S(t)  D(t) 


^       /   5dV  +  /   (u^P^^  +  Q^')dSj  =  y  £dV 
V<   D(t)      S(t)  D(t) 

Similar  integral  forms  hold  for  the  indii^idual  conserva- 
tion equations  (1). 
3.   Charge  and  Current 

We  note  the  conservation  of  charge 

^  'r  r 

and  define  the  current  densities 

J  =5'qv,J=yqu  =J  +  qu; 
u   ^  ^r  r'     ^  ^r  r    u   ^  ' 

J  is  the  current  observed  in  the  fixed  system  and  J 

"moving  with  the  fluid" J  the  convection  current  is  qu. 

Prom  (1),  the  individual  "conservation"  of  charge  is 

^X 
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whence,  adding, 

(9)  32  ^  ^i^  J  ^  0, 

Similarly,  from  (1)  we  can  obtain  individual  current 
equations 

(10)  Jr    (q  u^)  +  -K    (q  u^u^'  +  y  ^^  h    =  Y  (X^  +  (ft/). 
^   '      di     ^r  r    ^J^rrr    'rr     'rr    rr 

and  the  total  current  equation 

(11)  ^*  I  7^  ^%44  *  ^A'''   =  I  Yr"r  *  <^r4'- 

r  dx*^  r 

Prom  this  and 

(12)  ^  (qu  )  -.  u  -^  +  -^  (qu  uJ)  +  ^  _^  _  H  X 
follows     . 


^  ^  ^x«^    P  ^xJ    ^  ^r  ^  r    r  r    p 
The  current  equation, (11)  or  (13),  is  particularly 
significant  in  case  there  are  just  two  fluids,  one 
positive  and  one  negative.   In  this  case,  the  mass 
velocity  u  and  current  J  are  as  complete  a  description 
as  the  two  velocities  u-,  and  Up,  and  the  total  momentum 
equation  and  current  equation  can  be  used  to  replace  the 
two  momentum  equations.   Equation  (13)  can  be  written  in 

terms  of  J  and  u  alone  using  the  relation 
u  ^ 

(13a       T  q  vH''  =  -^^ /   ^    J^  JJ 

"^        ^  ^r  r  r   /      ^2       u  u 

(Y1-Y2)   P1P2 
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U»  Entropy 

A  local  the ririo dynamic  structure  can  be  introduced  by 
making  certain  assinnptions.   First  (this  has  already  been 
done),  extensive  quantities  such  as  mass,  momentum,  and 
energy  are  introduced  as  densities.   We  know  from  equili- 
brium thermodynamics  that  all  thermodynamic  parameters 
can  be  computed  from  the  single  function  "entropy  density" 
as  a  function  of  the  other  densities.   We  need  only  adopt 
the  same  fvinctional  relation  in  nonequilibrium  and  compute 
all  other  quantities,  e.g.,  temperat\are  and  pressiore,  by 
differentiation. 

We  define  the  entropy  per  unit  volume,  s,  by 
(li;)       Ts  =  e  +  ^  ^^p^   +  p 

and  write  the  second  law 
(15)       Tds  =  de  +  j;  tipdp^ 

and  the  Cxibbs-Duhem  relation 


0  =  ed(^)+Ip^d(^)  +  d(£). 


We  have  expressed  the  intensive  thermodynamic 
coordinates  per  unit  volume  rather  than  per  unit  mass 
since  this  provides  a  simpler  (additive)  transition  from 
individual  fluid  properties  to  total  fluid  properties. 
Expressed  in  terms  of  unit  mass,  (li|.)  and  (l5)  take  the 
form 

(16)     <   p   p   ^  ^^r  p    p 
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Prom  the  above  we  compute  that 

(17)  If  -^  div  (su  +  I  +  i  X  t^rPr^r)  =  ^  (£  -  u-X  +  Iix^cr  ) 

+  Q  .  grad(i)  +  Ip^v^  .  grad  (^)  -J  p^J  ^     . 

This  thermodynamic  formalism  is  appropriate  in  case 
it  is  sufficient  to  deal  with  the  total  momentvun  and 
energy  equations  alone,  although  taking  n  individual  mass 
equations.   In  case  we  wish  to  treat  all  n  velocities  and 
all  n  temperatures  as  primitive  variables,  it  is  necessary 
to  introduce  the  generalized  n  fluid  fon7ialism:-> 

r  r    r   ^r'^r    r 

(18)  "S   T  ds   =  de   +  u  dp 

'  r  r     r   "^r  ^r 

^  0  =  e  d(i-)  +  p  d(^)  +  d(^) 
r         r       r 

The  total  entropy  per  volume  is 

(19)  3^   =  I  »r 
We  compute  that 

(20)  ^^+  dlv  H   s^u^  ^l^)   =l^it^.   u^.X^) 

o- 

T 


+  Z  f^  (t^r  ""  i^r  ~  ""r*  '^r  ^  ■"  ^  ^  *  ^rad  (i-) 


V  1  niJ  ^^ 


■^     This  is  significant  only  for  perfect  gases,  i.e., 
when  there  is  no  significant  intercomponent  energy. 
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It  is  interesting  to  note  that  the  combination  u,   +  •:?  u 

^  *^r   2   r 

would  have  appeared  in  such  a  fashion  as  to  warrant  the 
name  chemical  potential  if  we  had  chosen  as  thermodynamic 
coordinates  mass,  total  energy,  and  momentum  instead  of 
mass  and  internal  energy,  (see  Grad,  J.  Phys.  Chem,,  56, 
1039-101+8,  1952). 

A  few  words  should  be  said  about  the  difference 
between  (1?)  and  (20).   First  of  all,  the  "entropy"  s-"- 
is  not  the  same  as  sj  specifically,  it  is  not  a  function 
of  the  same  variables.   In  the  one  fluid  formulation,  (17)» 
there  is  an  entropy  source  related  to  diffusion  velocities 
and  gradients  of  the  chemical  potentials,.  This  does  not 
appear  in  (20  )j  in  other  words,  the  individual  fluid 
velocities  are  here  treated  as  primitive  variables,  conse- 
quently reversible.   If  the  energy  sources  are  related  to 

the  forces  by  the  relation  o  =  u  •  X  ,  then  (20)  considers 
./  r    r    r 

this  as  reversible,  while  (17)  has  an  entropy  source  of 
magnitude 

g  -  u  .  X  =  5:  v^  .  X^  ; 
i.e.,  slipping  is  irreversible.   On  the  other  hand,  if  we 
have  energy  sources,  C'  »  in  the  absence  of  forces,  X  , 
and  have  £  =  ^^  -  0»  then  (1?)  considers  this  state  of 
affairs  as  reversible,  while  (20)  does  not  if  the  various 
fluids  are  at  different  temperatures. 

The  choice  of  the  entropy  function  s-"-  can  be 
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justified  on  the  basis  of  statistical  mechanics;  however, 
its  most  important  macroscopic  justification  is  that  we 
may  keep  the  concept  of  increasing  entropy  together  with 
a  fluid  description  involving  several  momentum  and  energy 
equations. 

As  an  example,  choose  for  the  force  X   and  energy  (> 
the  expressions 

r  X^  =  q^E  +  q^u^  x  B 

P^   =  u      '    X^   =   q  u^  .  E 
^  ut    r    r    r  r 

(E  and  B  are  the  electromagnetic  field  vectors, )   There 

is  no  entropy  production  from  this  source  in  the  second 

Case,  (20),  while  in  (17)  we  have 

(22)  C  -   u.X  =  J  .E   =  ( J  -  qu)  •  (E  +  uxB ) 

The  entropy  source  is  the  product  of  the  current  and  electric 
field  each  measured  in  a  system  moving  with  the  instantaneous 
local  fluid  velocity,  u. 

Another  possibility  is 

f  X^  =  q^E  +  (q^u^  +  J^)  X  B 

(23)  )     ^ 

(  C   =  E  .  (q  u   +  J  ) 
\  ^r       ^r  r    r 

The  extra  currents,  J  ,  could  be  introduced,  for  example, 

to  take  into  account  magnetic  polarization.   In  this  case, 

we  have 

^X  =  qE  +  JxB 

(24)  i    p 
O  s  E  •  J 
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by  redefining 

(Instead  of  just  ^  ^r^j,)*  whence 

£-  u  •  X  =  E   .J 

as  before,   but   also 

(26)  ^^  -  ^r   -   y^^  =  i^  +  n^  X  B)   •    J^     . 

It  will  be  observed  that  we  have  not  Included  electric 
or  magnetic  polarization  terms  in  the  thermodynamics. 
Nevertheless,  we  are  not  restricted  to  the  case  of  unit 
dielectric  constant  and  permeability.   To  achieve  this  we 
need  merely  adopt  the  viewpoint  ahat  J  and  q  include  polari- 
zation current  and  charge  respectively  (cf,  Stratton, 
equation  (7),  p.  12).   Moreover,  this  Interpretation  is 
almost  forced  upon  us  in  any  consideration  involving  energy. 
The  usual  formulation  in  which  electric  energy  is  represented 
by  jE'dD   rather  than  f^   E^  (Identical  remarks  hold  for  H 
and  B )  is  a  useful  expedient  only  when  the  internal  energy 
of  the  medium  is  disregarded.   In  certain  very  simple  situa- 
tions, it  is  convenient  to  distinguish  electric  charges 
which  are  \mder  the  direct  control  of  the  experimentor 
from  other  "polarization"  charges.   In  the  process  of  polari- 
zation, there  is  a  transfer  of  energy  between  the  electric 
field  and  the  dielectric  (Internal  thermal  energy).   If 
one  is  not  interested  in  this  transfer  of  energy  (which 
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is  frequently  reversible),  one  can  maintain  an  artificial 
conservation  of  those  energies  which  one  chooses  to 
acknowledge  by  using  /  E«dD  instead  of  / -i  E  •   In  oomplex 
problems  of  fluid  flow,  retention  of  D  and  H  would  seem 
to  be  a  hindrance  rather  than  a  help.   It  should  be 
remarked  that  equations  (8)  and  (10)  -  (13)  are  valid 
only  in  the  absence  of  polarization. 

In  succeeding  notes  chemical  terms  will  be  ignored 
unless  the  context  plainly  demands  their  inclusion. 
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